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I. 

INTRODUCTION 

In 1990, Ringwald and Espinosa |2j presented a calculation of possible instanton 
induced baryon number violation accompanied by enhanced production of a large number 
of Higgs and W,Z bosons at high energy^ ~ n m/g, where n is the number of outgoing 
particles of mass m, and g is a coupling constant of the theory). This has called attention 
to a more general and profound problem of qantum field theory. The problem is: for a given 
renormalized quantum field theory with weak renormalized coupling constant, what is the 
nature of the multiparticle production in that theory. 

It has long been known that the perturbation-theory expansion in a theory with weak 
coupling can fail in high orders because of the factorial growth of the coefficients in the 
series [[J. So this problem can not be solved perturbatively because the perturbation series 
diverges rapidly at high order. In fact, independent of the instanton calculations, Goldberg 
H, Cornwall|J and, later Voloshin|J, Argyres, Kleiss, and Papadopoulos 0,and Brown|8| 
pointed out that in 4 scalar field theory the contributions of just the tree graphs to the 
multiparticle production amplitudes gives the k\g k behavior. 

These complete tree graph calculations lead to a dramatic exponential growth of the 
rates for few — > many processes as the total energy E and number of final particles n 
increase together, and E exceeds some critical value of order E > nm/g. Something must 
turn this exponential growth off as it saturates unitarity bounds, but the phenomenological 
consequences could still be spectacular. Or perhaps the tree graph calculation is simply 
misleading - the dramatic exponential turn on only occurs for values of n ~ 1/g for which 
the perturbation series diverges. Perhaps the answer is given just by taking the first few 
(small) terms of an asymptotic perturbation series with a small coupling constant. On 
the other extreme, A number of semiclassical calculations have suggested that these high 
multiplicity events are exponentially suppressed. || 

It was found by Mawhinney and Willey |10[ that a nonperturbative approach could be 
made through Monte Carlo simulation on a lattice. It was shown that using the analytic 
continuation from Euclidean lattice theory to Minkowski space quantum field theory, one 
could turn spectral function sum rules into an upper bound on inclusive multiparticle pro- 
duction. The program was carried through for 4 theory in 1 + 1 dimension because, as 
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explained below, large lattices are required, and we only had the computing power for large 
two-dimensional (256 2 ) lattices. No signal for significant highmultiplicity processes was seen. 
However, the 1 + 1 dimensional theory is a very special case in field theory, (no angular mo- 
mentum, no concentration of radial wave function at spatial origin,..), so it is of interest to 
know the results for the 2 + 1 dimensional theory. The tree diagram considerations which 
suggest large 1 — > n amplitudes in 4 are mainly combinatoric, and hold as well in three 
dimensions, as in two (or four). Furthermore, the 4 theory is superrenormalizable in three 
dimensions as well as in two. This is an essential element of the analysis presented below. 
(In four dimensions one has to deal with the presumed triviality of the theory and a cutoff. 
In less than four dimensions, the existence of the theory as the analytic continuation of the 
continuum limit of the Euclidean lattice theory is established. [fLl) , [T2|| . ) 

To deal with large three dimensional lattices we built a Beowulf parallel supercomputer 
with sixteen CPUs. We also used a Message Passing Interface (MPI) function call to paral- 
lelize our simulation program. We are able to achieve overall 7.4 Gflops at High Performance 
Computing Linpack Benchmark(HPL). The theoretical peak performance for sixteen CPUs 
is 12 Gflops. Compared with a 16 Gflops CRAY C90 at 16 CPUs, our parallel cluster has 
much better price/performance ratio. (For details of the construction and performance of 
this Beowulf computer, [|l3"f ). 

II. 

SPECTRAL FUNCTION SUM RULES AND BOUND ON HIGH MULTIPLICITY 
PRODUCTION PROCESSES 

We start with some standard definitions to set the notation and display the ingredients 
used in the analysis below. In d spacetime dimensions, the canonical (unrenormalized) scalar 
field satisfies canonical commutation rules 



[0(z),0(O)] 



zo=0 



i8 d -\x) 



(1) 



The Wightman two-point function and its spectral decomposition are 




(2) 




(3) 
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2np(p 2 ) = J(d d x)e^(0\4>(x)m\0) (5) 
= (27r) d 5:/ d )(p- Pn )(O|0(O)|n)(n|0(O)|O) (6) 

The states \n) are members of a complete set of in-, or out- states, and are labeled as 
\q 1 , ...,q n ) with 

n 

Pn = E "(a) = ^ + ™ 2 (7) 

0=1 

The boldface sum over n implicitly includes some relativistic kinematical factors relating it 
to the relativistic invariant phase space integration. 

5»w = ioxoi + g - j n (/ i* - • ■ ■ < 8 > 

Relativistic invariant phase space is simply related to this 

W^-^S(/^)=/*.( £ ) («) 

where is the total cm energy E = YZ=i Q(a) 

Separate out the single particle contribution to the spectral function. 

p{p 2 ) = Z6(p 2 - m 2 ) + p{p 2 ) (10) 

and 



( 9 |0(o)|o) = Vz (ii) 

The CCR (HD imply the sum rule 

"dn 2 p(K 2 ) (12) 



The integral is convergent in less than four space-time dimensions. Together with the posi- 
tivity of the spectral function @, these equations imply 

< Z < 1 (13) 

For free field, Z — 1, and the field connects vacuum to single particle state only. 
The renormalized field (x) is defined by the normalization condition 

(gk/(0)|o> = l (14) 



Thus (j)(x) = \Z~Z(tf(x). 

For the multiparticle states (|)|),it is an exercise with the LSZ reduction formulas to obtain 

(ft . . . g„|0(O)|O) = J^G F (p 2 n )Z n r n+1 ( qi ...q n ) (15) 

Here r n+ i(gi . . . q n ) is the Fourier transform of the vacuum expectation value of the T- 
product of n+1 <ft fields, with the energy- momentum conserving delta function left off. The 
underline indicates that the external lines are all amputated (qi . . . q n are all on-shell.) Gf 
is the complete unrenormalized time-ordered two-point function. 
The renormalization of all these quantities is 

P = Zrf, g f = ZG F , r (n+1) = (-^±r) V(n+i) ( 16 ) 

Then (||) becomes 

27r// w = / d$ n (E)\G- F (pl)\ 2 \i {n+i) ( qi . . .g n )| 2 (17) 

The G F depends only on the total four-momentum p n (E in the c.m.)which is fixed in the n 
particle phase space integral, so it can be divided out. Divide also by particle density factor 
2E 

I^H = IE /*>WIV,„ +11 ta - «.)!' = r w W (is) 

which is the decay rate for one variable mass (E) <ft particle to decay into n on-shell (f) 
particles. The totally inclusive rate is 

r(£) = ]Tr (n) (£) (19) 

n 

One can think of the single timelike 4> particle as formed by the annihilation of an // 
pair, coupled weakly to the <fi field. Then up to trivial kinematic factors, this is the analog 
of the QCD ratio R of hadron production to n + fi~ production in e + e~ annihlation. We do 
not write a slash for T. There is no reason to introduce a bare rate, so we simply write V 
with no slash for the physical quantity. 

We now relate this rate to the inverse two-point function. (This has a simple diagra- 
matic representation. The probability for one timelike <j) particle to go to many on-shell <fi 
particles is given by the Cutkosky rules applied to the cuts of the self-energy diagrams). 
The spectral representation of the Feynman function follows from (2) 

Gf( P 2 ) = ~^— 2 + / ^ 2 2 ^ K l\ (2°) 

pz _ m z j pz _ K z _|_ le 
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The imaginary part is 

%G F (p 2 ) = -Kf>{ P 2 ) (21) 

The imaginary part of Gp 1 is 

= = (22 > 



From all the information about Gp, we can write a dispersion integral for Gp 1 |T4] 



with 



7(^ 2 ) 

^(p 2 ) = vr 7 (p 2 ) = |r(£), (p 2 = £ 2 ) (24) 



We note that all of the properties of Gp enumerated so far are consistent with the possible 
existence of zeros of Gp on the real axis. If Gf has no additional poles, only one zero is 
possible, on the real axis in the gap between the pole and the start of the continuum. In 
any case, one can show that if such zero(s) of Gp exist, positivity of the spectral function is 
sufficient to show that the bound derived below is only strengthened [[L3|] 
The analytic continuation of the two-point function is immediate 

P° -> m, P 2 M = -pjs, g f(Pm) = -G E (p%) (25) 

GiV ) = (p* + m")(l + / Jig,^ ) <*) 

In this and subsequent equations in this section p 2 is Euclidean p\ The desired sum rules 
are now obtained by expanding both sides of this equation in powers of p 2 and equating the 
coeficients. The expansion of the left hand side defines the quantities Z ,m 

G E \p 2 ) = ^(m' 2 +p 2 + ...) (27) 

while the right hand side is 

G E \p 2 ) = m 2 (l + / dK 2 , 2 1{K \ 2 ) +p 2 (l + / d. 2 ^) + . . . (28) 
J [K z — m z )K z J K 

Matching coeficientss of p 2 gives 
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Finally 

If the probability for production of high multiplicity states of cm. energy E turns on 
exponentially at some critical energy E+, the integral of T(E)/E 2 will turn on exponentially 
as the range of integration over E extends beyond E+. But by above, this integral is bounded 
by — 1, which can be extrapolated from lattice MC calculation of the slope of the inverse 
of the Fourier Transform of the two-point correlation. 



III. 

CRITICAL ENERGY AND LATTICE MOMENTA 



According to the papers cited in the introduction, there may exist a critical energy E+ 
(~ n m/g) at which the probability for production of high multiplicity final states suddenly 
begins exponential growth. We require an estimate of that critical energy and its relation 
to the momentum of the lattice propagator. Repeat the spectral representation for Ge(p 2 ) 
(p 2 is the Euclidean momentum squared). 

G ^=^h?+J d « 2 M (3i) 

k 2 is the square of the cm. energy of the state |n) contributing to the spectral function p 
(section two). For large p 2 ,the (convergent) integral is dominated by k 2 < p 2 . Then for the 
integral to be sensitive to dramatic behavior for k 2 « E 2 requires p 2 > E 2 . But the relevent 
p is the Euclidean lattice momentum which lies in the Brillouin zone — 7r/2 < p^ < n/2, or 
(lattice) p 2 < d(n /a) 2 . Combining these observations leads to the desideratum 



7T 



a 



2 



d - >Ei (32) 



Take the square root and divide by mass m 

E*/m < Vdn^ L (33) 

where I /ma = is the correlation length in lattice units. 

The estimation of for the tree diagrams of 4 in three plus one dimensions has been 
refined in the papers cited in the introduction. We have adapted these calculations to two 



plus one dimensions and find [13] 



E± m , „ 

-± = 74tt— 34 
m A 



Substituting this result into (|33|) we obtain 

74vr^ < y/dn£ L (35) 
A 

There is also a pair of general conditions for extrapolating results obtained from a cal- 
culation on a lattice to results in continuum limit: The correlation length should be much 
greater than the lattice spacing and much less than the linear lattice size. Thus 

74 TYl 

1 <<: TIT <h« N ( 36 ) 

Note that for fixed — of order one, these conditions can all be satisfied by a large enough 
lattice (N of order one or two hundred). 



IV. 

LATTICE SIMULATION AND CONTINUUM LIMIT 

The Monte Carlo simulation for the scalar field on the lattice is standard. The physical 
field, masses, and coupling constant are dimensionful quantities. When put on a lattice the 
lattice spacing a can be scaled out and the dimensionless action written in terms of dimen- 
sionless lattice fields and mass and coupling constant. The distinction becomes important 
in taking the continuum limit. In three (Euclidean) dimensions the scalings are 

= ~^=<t>L, (4 = \hI l , , A = ~X L (37) 

The lattice action is 




We start with a standard Metropolis update algorithm, and with each sweep of the lattice 
measure the expectation of </> and the two-point correlation. (For this discussion of the lattice 
calculation all quantities are in lattice units. We supress the subscript L until we come back 
to the continuum limit). We take the lattice Fourier transform of <fi and compute the output 
lattice propagator 

G(k) = ((/)(k)(f)(—k)) (39) 
and plot the inverse of this against the inverse of the free lattice propagator 

G \k) = k 2 + m 2 (40) 
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with 

P = 4^sin 2 (V2) (41) 

This is the inverse of the massless free lattice propagator. 

Z' and m' are obtained from the plot of G~ l (k 2 ) against the inverse of the free massless 
lattice propagator. (p7|) 

As we approach the contiuum limit we are also approaching the critical line fig (|J), and 
the simple Metropolis update begins to suffer critical slowing down. As one approaches the 
critical line, both the correlation length and autocorrelation time diverge. The relation of 



the correlation length £ and autocorrelation time r isJTR 121 



r ~ £ z ~ L d (42) 

where z is the dynamic exponent and r is measured in number of sweeps. In a system of 
finite size, the correlation length never really diverges. Cutting it off at the size of the lattice 
, the autocorrelation time increases as the lattice size increases as 

r ~ L* (43) 

when close to critical. Depending on the value of z, one can require a large amount of CPU 
time to generate a fixed number of independent data as one approaches critical on large 
lattices. For straight Metropolis update, z is close to two. See Table I. 

To overcome the critical slowing down, one must define appropriate nonlocal (or collec- 
tive) variables and a new dynamics for driving them. Progress has been made with some 
nonlocal algorithms for both discrete spin models and continuum fields. Swendsen and 
Wang [[16|j have used the Fortuin-Kastelyn |lj]] percolation map for the Potts model to define 



collective coordinates that allow domains to be inverted with zero free energy cost. Some 
modifications were proposed by Wolff fT%| 



The algorithm for the 4 field theory is a little more complicated. Brower and Tamayo|20|] 
proposed an algorithm for 4 field theory based on the Swendsen and Wang algorithm. Here 
we will follow Brower and Tamayo's method but switch to the Wolff algorithm. Table | com- 
pares the dynamic exponent values for the Ising model in different dimensions with different 
algorithms. One can see that the Wolff algorithm is more effective than the Swendsen and 
Wang algorithm in three dimension. 
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TABLE I: Comparison of the values of the dynamic exponent z for different algorithms for Ising 



model in various dimensions! 22 1. 



dimension d 


Metropolis 


Wolff 


Swendsen-Wang 


2 


2.167 ± 0.001 


0.25 ± 0.01 


0.25 ± 0.01 


3 


2.02 ± 0.02 


0.33 ± 0.01 


0.54 ± 0.02 


4 




0.25 ± 0.01 


0.86 ± 0.02 



In our simulation, the update algorithm consists of two parts: (1) a conventional Metropo- 
lis Monte Carlo update for the <f>(x) field and (2) Wolff cluster identification and flipping. 
Identification of a cluster of field, requires introduction of a discrete variable 

0(n) = s H |0(n)| (44) 

where = ±1. 

For details of the cluster algorithm see [13j Here we only present some illustration of its 
effectivness and correctness. In fig(|l] ) we show plots of autocorrelation in measurements of 
(0) for one set of input parameters, with and without the cluster algorithm. We see an 
improvement by a factor of roughly 200 in the number of sweeps required to reach e~l fall 
off. 

There is a price to be paid for this improvement in generation of uncorrelated equilibrium 
configurations. That is in the actual measurement of (0) itself. With only Metropolis update, 
on large lattices, there are very many sweeps between tunneling, and measurement of ((f)) 
is trivial. The cluster update changes sign of (0) frequently and it is harder to dig out its 
infinite volume (no tunneling) limit. 
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FIG. 1: The autocorrelation for the (0) with and without the cluster update algorithm. The e -1 
point for the one without the cluster update algorithm is at about 400(top). The autocorrelation 
length for the one with cluster update algorithm is about 2 only (bottom). 
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FIG. 2: Result of Monte Carlo evaluation of G~ 1 (k) for A = 0,/iQ = 0.001 plotted against exact 
analytic massless propagator k 2 = 4 J2 V sin 2 (k u /2) on a 32 3 lattice. Error bars on Monte Carlo 
results are too small to show up on this plot. 



An important check is to do the MC simulation for the coupling constant A equal to zero. 
For this case we know the exact analytic (free field) solution. In particular, mf" is equal to 
input bare mass squared (necessarily positive) and Z' = Z = 1) . This is shown in fig^ where 
G^ 1 from the simulation (with A = 0)is ploted against the free massless inverse propagator 
(A; 2 ). The result is a straight line of unit slope, offset by the mass squared (0.001 and not 
visible in this example). 

Although the physics is entirely dependent on A ^ 0, for the MC simulation program and 
the program for analysis of the MC data, zero is just one possible value of A. So the above 
result does check a substantial part of the simulation and data analysis programs. 



12 



For A 7^ 0, we have also, for a few sets of input parameters, paid the price in longer CPU 
time required when the cluster algorithm is turned off, and checked that the results (m f , 1/Z') 
are the same, within small statistical errors, with or without the cluster algorithm. 

We will return to another test, the Ising limit, after our discussion of the phase structure 
and the continuum limit. 

There are both general theoretical and heueristic and practical considerations involved 
in extracting real world results from the Euclidean lattice results. The general theoretical 
underpinning comes from the program of constructive field theory which has provided the 
results that for d < 4, the <fi\ Quantum Field Theory can be obtained as the analytic con- 
tinuation of the Euclidean 0^ field theory, which in turn may be obtained as the continuum 
limit of a lattice theory. This program is described in detail in two long and mathematical 
monographs ||Tl|,||12||. A simple example is again provided by the free theory (A = 0) which 



is completely determined by the two-point function . When the lattice propagator ([11]) is 
rewritten in terms of physical dimensionful quantities with the lattice spacing a inserted 
where required to give the propagator its correct physical dimension, the limit a — > is 
simply seen to be the free Euclidean propagator. However, in the general case, the output 
from the MC simulation is numerical, and the limit has to be done by extrapolation. To see 
what is involved, we consider the phase diagram fig (|3|). 

The 4 theory on a lattice of lattice spacing a has two dimensionless parameters, /Zq l , A^ 
(|57|). Since A is finite and /Zq is only linearly divergent (~ 1/a), both of these dimensionless 
parameters go to zero in the limit a — > O.Thus the point at the origin in fig(§) corresponds 
to the Euclidean field theory. So start on dimensionless lattice (iV 3 ) at some point in /Xq , A^ 
plane. Choice of parameters is subject to requirement that output £' L = l/m' L satisfy 

1 « & < N (45) 

This is an obvious condition for lattice spacing and volume effects to be small. To proceed 
toward the origin, go to lattice N 2 = 2N, and adjust input /!q l ,Al to get output £' L = 
2£' L . This corresponds to keeping same physical volume and correlation length while reducing 
lattice spacing to half of starting value. A sequence of these steps corresponds to one of 
the dashed lines in fig(^). Note that at each step £' L is increasing (doubling) so the first 
inequality in ( [45] ) is increasingly better satisfied (The second remains the same). 

The 4 field theory (a — ► 0) is specified by a single finite dimensionless constant, which we 



13 



Ai 




FIG. 3: The phase diagram for the three-dimensional, single component scalar </> 4 theory. The 
solid line separates the symmetric phase from the broken symmetry phase. The dashed lines are 
lines of constant physics. 

may take to be X/m'; so "constant physics" is implemented by fixed value of A/m'.Also, the 
two other finite dimensionless quantities we can measure, 1/Z', and ((f)) / y/rri/ ,are functions 
of X/m' in the limit of a — > 0. As the origin is approached along one of the dashed lines, the 
values of these quantities should approach their continuum values. In practice, When the 
change in values in succesive steps is less than the estimated error, we accept thcie values. 

We briefly discuss the relation of this work to the Ising model of Statistical Mechanics. 
This is in itself of some interest, and it also provides another check of our MC simulations. 
The EFT limit is reached by following one of the dashed lines into the origin (Al — * 0). The 
Ising model limit is approached by following a dashed line (in broken symmetry phase) out 
to infinity (Al — > oo).The Euclidean lattice Lagrangian Q38D can be rewritten as 
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TABLE II: The Monte Carlo simulation of 32 3 ^3 theory at Ising limit(A — > 00). The table shows 
the critical temperature (3 C at different A values. 



A 




A 




100 


-31 


.31 


.25 


300 


-75 


.25 


.23 


1000 


-228 


.228 


.222 



c = - E m0) + E \ (f® - -^f^J ( 46 ) 

All quantities here are in lattice units (subscript L suppresed).The sum over angle bracket 
is standard notation for sum over nearest neighbors, each pair counted once. An irelevent 
constant has been added. Compare with the Hamiltonian of Ising model 

(13) 



The limit A — > 00 will exponentially suppress the second term in equation fl46|) , and convert 
the continuous field variable to a spin variable. 



with <7? = il and 



lim^(i)=±i/-^±^ = ^ (48) 

A^oo V A v 



^ = ~^ + 2d (49) 
A 



To compare with the Ising limit, we do the simulation with the A ^> 1 and calculate the 
VEV (0) for different /j,q. We start in the symmetric phase with one A value and vary the 
value of /ip until the VEV is no longer zero. In our simulation, when A ^> 1, we get a very 
sharp transition from the symmetric phase to broken symmetry phase. So, it is very easy 
to identify the critical $ . Table |I| shows the (3 C from the simulation at 32 3 with different 
A. The third row shows that when A = 1000, the j3 c from our simulation with a continuous 
field variable is very close to the (numerically) known value 0.2217... for the 3d Ising Model 
This give us additional confidence in our MC simulation. 
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V. 

ANALYSIS AND RESULTS 



A fortran program mplementing the Metropolis and cluster update algorithms described 
above provides the equilibrium lattice configurations i.e. the sets of values 4>{n) t for t from 
1 to Nt where t denotes the t'th sweep through the lattice andiVt is the number of sweeps 
through the lattice after thermalization. For large lattices with Nt less than the time before 
tunneling, and the cluster update turned off, the vev is 

Nt 

<<P>=l/N t l/N^H^)t (50) 

1 n 

When the cluster update is turned on, the sign of the vev as defined above changes 
frquently and < > is very small, even deep in the broken symmetry phase. So the vev is 
then defined as 

N t 

<|0|>=l/iV t l/iV 3 ^|^0(n) t | (51) 

l ft 

The program also does the lattice Fourier transform of and constructs the k-space 
lattice propagator. This output is then input for a Mathematica program which determines 
m' 2 /Z' and l/Z' from fit of output lattice propagator 

G(P)- 1 = (l/Z') (m /2 + P + . . .) (52) 

The estimated errors asociated with l/Z' and ml for fixed lattice (N) are statistical, origi- 
nating from the Monte Carlo, and systematic from extrapolation to k equal to zero in the 
fit. (We can not use the data point k = because it contains a volume singularity) The 
results are collected in Table [TTT] 

We can see in the Table the implementation of the approach to the continuum limit by 
a sequence of steps in which the physical correlation length (physical mass) is kept fixed 
while the linear number of lattice points (N) and the correlation length in lattice units (£' L ) 
are both doubled. The ratio N/£' L (£/£) is kept fixed (~ 5.6). We also see that X/m' 
keeps the same value within expected errors through these steps as do the dimensionless 
quantities, l/Z', (0) / \fm , consistent with their being functions of X/m' only. This suggests 
that these values are insensitive to lattice artifacts and approaching the continuum limit. 
This is illustrated in figures (§J),(f|), (f|). 
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TABLE III: The results of fit to data from the simulations on different size lattices(iV 3 ),the in- 
put parameters //q , A£,and the output dimensionless coupling parameter A/W,and the output 
correlation length £' L = l/m' L ,the slope 1/Z', and the dimensionless ratio < \c/)\ > j\fm! . 



N 3 




Al 


A 

m' 


e L 


l 

Z> 


M 

\fm! 


32 3 


-0.161 


0.2 


1.13(2) 


5.7(1) 


1.008(2) 


0.78(1) 


64 3 


-0.078 


0.1 


1.13(4) 


11.3(4) 


1.012(3) 


0.80(2) 


128 3 


-0.0384 


0.05 


1.11(5) 


22.11(1.0) 


1.012(5) 


0.80(2) 


256 3 


-0.01905 


0.025 


1.13(8) 


45.4(3.0) 


1.010(8) 


0.79(3) 



We also check the absence of finite volumn effects. We hold the input 
parameters, [1q , A^, fixed and vary N. In Table IV we see that when the linear size of the 



lattice is greater than three times the output correlation length, the output quantities are 
unchanged under further increase in N. Note that in Table |IT| the ratio of N to £' L is kept 
fixed at 5.6 as one moves along a dashed curve in fig(0) 

TABLE IV: The finite size effect with fixed input parameters $ = 0.01905, \ L = 0.025. 





Al 


N 3 




m 1 


l 

Z 1 


(M 

\J m' 


-0.01905 


0.025 


64 3 


41.5(6.5) 


1.04(16) 


1.002(3) 


0.78(14) 


-0.01905 


0.025 


128 3 


45.1(3.8) 


1.127(97) 


1.0011(35) 


0.77(10) 


-0.01905 


0.025 


256 3 


45.4(3.0) 


1.13(8) 


1.010(8) 


0.795(97) 
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Continuum Limit 

lfZ' 



1.05 



□ lfZ' 



[] [] 



§ 1 - 



t T T | 

N=32 N=64 N=128 N=256 

0.95 ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

100 200 300 

N 



FIG. 4: Continuum limit for yj. The plot shows the insensitivity to lattice spacing for yj. 
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FIG. 5: Continuum limit for the A. The plot shows the insensitivity to lattice spacing for A 
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FIG. 6: Continuum limit for the H=. The plot shows the insensitivity to lattice spacing for 
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The significance of these results is brought out in the qualitative behavior of the quantities 
in Table |V[ 

TABLE V: Size of lattice and condition for appearance of high multiplicity inelastic processes 



1 « 


74 m' 


Cl < 


< N 


l/Z' - 1 


1 


37.8 


5.7 


32 


.008(2) 


1 


37.8 


11.3 


64 


.012(3) 


1 


38.5 


22.1 


128 


.012(5) 


1 


37.8 


45.4 


256 


.010(8) 



The condition (|36| 



74 m' 

1<< 7ST <(L<<N 



(53) 



is not satisfied for the first three lattices, but is clearly satisfied for the last (finest grained) 
lattice. And there is no significant change in the value of l/Z' — 1 in going from the first 
lattice to the last. We consider this to be strong nonperturbative evidence that there is no 
exponential turn on of high multiplicity production processes for energies greater than some 
critical energy of order 

Allthough the most significant feature of TablejV] is the absence of large change in l/Z' — 1 
as one passes into a region of parameter space in which the divergence of perturbation theory 
suggests the possibility of exponential growth of high multiplicity processes; one can still 
ask: what would be a large value or a small value of l/Z' — 1 in this context? What is the 
significance of the single number 0.010 ? In the weak coupling regime, we can take low order 
perturbative effects to define small. We have calculated l/Z' — 1 through two-loop order, in 
broken symmetry phase. ||13|| . We have chosen to specify the constructed cj)\ theory by the 

^-7 The result is 

1_ . A A 
Z~' 



output A/to'. For the perturbation expansion parameter (g) we take g __ _ 



1 + 



-7(3/4) + 



! (-1.83) 



(54) 



OTrm' 87rm' 

For our particular lattice MC solution, A/to' = 1.13 (g = .045). This gives the perturbative 
value 

l/Z' - 1 = . 0337 -.0037 + ...-. 030 (55) 
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The lattice MC result from Table [TV] is .010(8), which is close to consistent with the per- 
turbative estimate. So any nonperturbative contribution is of the same order as or smaller 
than the perturbative estimate, which is small. 

We can still ask: If there were important high multiplicity processes contributing to (ftO"D, 
how would they show up in this analysis? From all the spectral equations in section II, we 
see that this would show up as Z substantially less than one. In this case, the last step in 
(j30|) leads to an inefficient bound. So we go back to 

2/ dB £W = Z(l-l)<l-l (56) 
7r J E 2 Z Z y 1 

By manipulation of the equations of section II, one can show that ZjZ' < 1. Then on the 
righthand side of the equation 

Z(i - 1) = J? - Z < 1 - Z < 1 (57) 

So even if the decay were dominated by high multiplicity final states, the integral of the 
inclusive rate is bounded by one. 

(58) 

While in this limit we expect 1/Z' to be much greater than one i.e. 1 /Z' -1 is a very inefficient 
bound in this case. 

It is also interesting to compare the perturbative and nonperturbative calculations of 
< |0| > /^/m/.The perturbative expression is 

For X/m' = 1.13 this is 0.81. The lattice MC value from Table-TV] is 0.795(95) ( or rounding 
0.80(1) ) 

In conclusion: The absence of any significant change in the value of 1/Z' — 1 from the 
lattice MC calculation as one moves from coarse grained to fine grained lattices, and that 
the value is close to the two- loop perturbative result, are strong evidence that there is no 
exponential enhancement of the rate for the high multiplicity processes 1 — > n for n ~ 1/g. 
We are not sensitive to the presence or absence of exponential suppression of rates for 
asymptotic values of E,n, which by definition are very small and only effect terms far out in 



the series fl54|) 
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